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Abstract. 

We check how the change in surface conditions between the Sun and red giant branch stars changes the 
characteristic surface convection length scale to be used in models. We investigate the question in the case 
of the mixing length theory and of the phenomenology of full spectrum of turbulence. For the observational 
part, we rely on independent measurements of effective temperatures and interferometric radii of nearby red 
giants. We find that the local red giant branch cannot be explained taking into account the solar calibrated 
convective length scale. 
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1 Introduction 

In low mass stars the effective temperature and the radius estimate the efficiency of surface convection. The 
Sun and the red giant branch stars (hereafter RGBs) have very different surface conditions. A red giant surface 
gravity and energy flux are much smaller than the solar ones. The purpose of the work we present is to check 
how this changes the surface convection efficiency. First we build calibrated solar models using the mixing 
length theory (hereafter MLT) or the phenomenology of full spectrum of turbulence by Canuto, Goldman & 
Mazzitelli (1996) (hereafter CGM) for surface convection. Then relying on the same input physics, we build 
red giant models. The red giant models radii and effective temperatures are compared to an observational 
sample of 38 objects for which the radii are known directly through interferometry to better than 10 percent. 
Absolute luminosities and effective temperatures of these objects are also accurately known. In the next section 
we give the main ingredients of the models affecting the radius and effective temperature and address the solar 
calibration. In section|2]we describe the RGBs giant sample. Section 2] investigates the changes in characteristic 
convection length scale from the Sun to RGBs. We conclude in section [5] 

2 Models inputs and solar convection length scale 

The radii of solar and RGBs models are tuned by the outer thermal gradient. It depends on four ing redients: 



3 radn oi solar and KLrOS models are tuned by the outer tnermai gradient, it depends on tour ing redients: 
i) The opacities. We rely on the OPAL opacities and, below 5600K, on the Ferguson et al. (l2005l) opacities 



The composition is either assumed to be solar with X=0.7392, Z= 0.0122 (Asplund et al. ( 20051 )) or slightly 



subsolar at [Fe/H]=-0.17. The metal repartition is always the one of lAsplund et al. I (|2005l ). 

ii) Convection efficiency to transport energy. The poor efficiency of outer convection induces the subsurface 
superadiabatic gradient which in turn sets the entropy of the deep convection zone. A lower entropy level of the 
deep convection zone means a less dense envelope and thus a wider radius and a lower effective temperature. 
To mo del the inefficien t conv ection we use two simplified local treat ments of convection: the mixing lengt h 



theory ( Bohm- Vitense (1958)) and the full spectrum of turbulence of Canuto, Goldman fc Mazzitelli ( 1996f) 



For both treatments the characteristic convection length scale A is assumed to be a constant fraction of the 
local pressure scale height: A = aHp. 
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iii) Atmospheric effects. At the edge of the star the diffusion approximation does not hold for photons 
and lines strongly affect the radiative transfer: these effects are expressed by the temperature-optical depth 
relations that are provided by atmosphere models. We use two series of non-grey atmosphere models as outer 
boundary conditions. The first series of relations (T(t)^ = T|Jgfgrid(T)) is computed with the PHOENIX/ID 
atmosphere code where the convection is handled using the M LT. The second series of temperature-optical 
depth relations is computed with the Atlasl2 atmosphere code ( Castellil ( 20051 )). We modified Atlasl2 to use 



the CGM prescription. Each type of atmosphere models are used with the corresp onding phen omenology of 
convection in the deeper regions as is necessary for consistency in the models (Monta lban et alTI ([2004)). 

iv) The equation of state influences the radius through the adiabatic exponents. We use the OPAL EoS. 

We assume that L0 = 3.846 lO^'^erg.s"^ and Kq — 6.9599 lO^'^cm and begin the solar evolution on the zero 
age main sequence. The calibration in luminosity, radius, and metal-to-hydrogcn ratio are achieved to 

better than 10^"* at the age of 4.6 Gyr for both MLT and CGM convection prescriptions. In the MLT framework 
we obtain Q;mit0=l-98, in the CGM framework we obtain Q;cgm0=O.77. 

3 The red giant sample 



We first queried the CHARM2 catalogue (jRichichi et al. I (|2005l )) to obtain all direct measurements of giant 



and subgiant angular diameters up to 2004, with effective temperatures in the range from 5000 K to 5500 K. 
We then searched the literature for more recent observations, and added the measurement of 7 Sge, S Eri, 
^ Hya, and the recent high accuracy CHARA/FLUOR measurements of eOph and rySer. The conversion of 
uniform disk angula r diameters to limb-darkened values was done usin g linear l i mb-d arkening coefHcients by 
IClaret et al. ] (Il995h . which are based on stellar atmosphere models by iKurucz I (jl993l ). Our sample contains 



38 giant and subgiant stars with spectral types from G5 to MO. The distances to the selected stars range from 
11 to 110 pc. Thanks to this proximity, we neglected the interstellar reddening for the computation of the 
bolometric luminosity. The accurate parallaxes and interfcrometric angular radii estimates allow the objects to 
have a relative uncertainty in the linear radius smaller than 10%. The average metallicity is slightly subsolar 
[Fe/H] = —0.17 with no object below -0.44 and no object above 0.04 but one exception at 0.13. 

4 Red giant branch calibration 

We model RGB stars up to 10'^ solar luminosity with exactly the same physics as in the solar mode l s. Th e 
microscopic diffusion is accounted for in any model warmer than 5000K following IP^ofRtt fc Michaud I (|l993[) . 



This is important in order to obtain correct ages as diffusion speeds up the main sequence evolution. After the 
first dredge-up though (Toff < 5000K) diffusion effects become negligible. Provided the atmosphere boundary 
models and the opacity tables are unchanged, there are three main models inputs that change the position of the 
RGB: the mass, the metallicity and the surface convection characteristic length scale. The latter parameter is 
what we want to constrain. As mentioned above the more efficient the convection, the smaller the radius and the 
higher the effective temperature at a given luminosity. The metallicity of the sample is known and is therefore 
no hurdle. The masses of the stars however are not known: unlike RGB stars of a globular or a Galactic cluster, 
the RGB stars of the sample are field stars that presumably have different masses and ages. Yet it is possible 
to set an upper limit to the age of local red giants or equivalently a lower limit to their masses. The limit is 
given by the age of the Galactic disk and the evolutionary ti mescale of its low mass stars: for objects in the 
slightly subsolar metallicity range (-0.25 < [Fe/H] < -0.14). ILiu fc Chabover I (I2OOOI) suggest a maximum age 



of 11.7 ± 1.9 Gyr. In this study, we consider models that have reached lO^L© on the RGB by that age as our 
RGB stars exhibits [Fe/H] = —0.17 on average. We can s et broader upper l imits to the local Galactic disk age: 
it is certainly younger than the Universe 13.7 ± 0.1 Gyr ( Komatsu et al. ( 2009f )). In the next subsections we 



focus on the lower envelope of the RGB i.e. on these stars with the lower effective temperature or larger radii 
at a given luminosity (see Figure [T]). They are also the oldest and lowest mass stars of the sample. 

4.1 Mixing length theory 

Figure [2] left panel shows the six lower envelope stars of the RGB and features different evolutionary tracks. 
Let's define x? = w[ ^Tt^^°°" ^^^^^ ^^s"^ ^'^'^ ^^"s" = 130K are respectively the observed 

Teff of an object, the Tos of the model having the same luminosity as the object and the uncertainty on the 
observed Teg. N is the number of objects considered. Figure [2] left panel models in solid line, dashed line and 



Solar and RGB surface convection 
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Fig. 1. HR diagram of the RGB stars sample. The triangles define the lower envelope of the RGB in the sense that they 
show the stars with the lowest effective temperature for a given luminosity. 



dotted line all use the solar calibrated value of amit = 1-98. The solid line track corresponds to a O.95M0 star 
reaching lO'^L© at wll.5 Gyr. This model has [Fe/H] = —0.17 and an helium fraction Y=0.2582. It is clearly 
too warm to fit the lower envelope of the RGB: Xi=3.8. The dashed and dotted lines correspond to less massive 
and helium poorer stars with O.9M0 and Y=0.2482 respectively. The lower mass model is extreme in the sense 
that it reaches lO'^L© at «13.9 Gyr (or lO^L© at ?al3.88 Gyr), which is older than the current age estimate of 
the Universe. The heli um poor model is also extreme in the sense that its helium fraction nearly is that of Big 



Bang nucleosynthesis (jCoc et al. I (|2004l )) and evidently cannot be lowered any further. Both two last models 
are in slightly better agreement with the data than the former one. Yet they do not provide a good fit to the 
observations. This demonstrates that mass or helium fraction cannot be changed to improve the agreement to 
the observations. Models with lower a,nit than the solar value provide much better fits to the data. The three 
dotted-dashed line is the track of the model with a^ait = 1-68, O.95M0, and [Fe/H] = —0.17. This model reaches 
IO^Lq at 11.5 Gyr and has x?=0.16. The long-dashed line model has a^it = 1-68, 1.13Mq, and [Fe/H] = 0. 
It reaches at 7.5 Gyr and has Xi=0.078. The analysis in the HR diagram suggests a smaller than solar 

calibrated characteristic length scale for the MLT. 



4.2 Full spectrum of turbulence phenomenology 

We will now perform a similar analysis as above but regarding the CGM phenomenology. Furthermore, instead 
of using Teff we will use the interferometric radii. We therefore do not set ourselves in the HR diagram but in a 
luminosity vs. square of radius diagram (see Figure fright panel). Let's define xl = Ei=i w[ ar^ ^^^^ 
again N is the number of objects considered. Robsi R-mod and ARobs are respectively the observed radius of an 
object, the radius of the model having the same luminosity as the object and the uncertainty on the observed 
radius. Figure [2] right panel shows the nine stars with the largest radii at given luminosities. They were selected 
by considering the nine largest deviations to a linear fit of the whole sample in the luminosity vs. square radius 
diagram. The solid line track in Figure [2] right panel corresponds to a O.95M0 star reaching lO^Lo at wll.6 
Gyr. This model has [Fe/H] = —0.17, an helium fraction Y=0.2582 and the solar-calibrated value acgm = 0.77. 
This model corresponds to too small radii to fit the lower envelope of the RGB: x| = 7.2. As in the case of 
the MLT, changes in mass and helium fraction are unable to significantly improve the agreement to the data 
and we will not discuss them. On the opposite if we decrease acgm down to 0.62 we recover a good agreement 
to the observations. The three dotted-dashed line is the track of the model with acgm = 0.62, O.95M0, and 
[Fe/H] = —0.17. This model reaches lO^L© at 11.8 Gyr and has X2=0-'''0. The long-dashed line model has 
a^it = 0.62, 1.17M0, and [Fe/H] = 0. It reaches IO^Lq at 6.9 Gyr and has X2=0-40. The analysis in the 
luminosity radius diagram suggests a smaller than solar calibrated characteristic length scale for the CGM. 



5 Conclusion 



We modelled the Sun and local RGB stars in order to check if the change in surface conditions implies a change 
of the characteristic convection length scale A for two local treatments of convection: the MLT and the CGM. 
In both cases we assume A = aHp. Therefore we do not consider the original version of the CGM where the 
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Fig. 2. Left: Position of the six lower envelope stars of the sample with Tcff errorbars and various evolutionary tracks 
computed using the MLT. See text for a description of the corresponding models. Right: Position of the nine stars of 
the sample with largest radii and errorbars. The various evolutionary tracks in overplot were computed using the CGM. 
See text for a description of the corresponding models. 



characteristic convection length scale is the distance to the boundary with the region stable with respect to 
convection. At a given absolute luminosity a tunes the position of the RGB in effective temperature and in 
radius. We have accurate data on absolute luminosities, effective temperatures and radii of all RGB stars we 
consider. The location of the coolest stars or equivalently the largest radii stars of the sample suggest a decrease 
in surface characteristic length scale with respect to its solar calibrated value. We have shown the decrease to 
be required for the MLT in the HR diagram and for the CGM in the luminosity radius diagram. However we 
could have inverted the diagrams with res pect to the c o nvect ion treatments, the result would have been similar. 
The reader will find many more details in iPiau et al. I (|201l[ ) where we also specifically address the three RGB 
stars of the sample with asteroseismic mass estimates. The combination of interferometric and asteroseimic 
data clearly opens up new perspectives in the understanding of stellar fundamental parameters and how they 
can be used to constrain stellar physics (jHuber et al. I (|201l[ )). 
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